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Abstract 

The production of Hawking radiation by a single horizon is not depen- 
dent on the high-frequency dispersion relation of the radiated field. When 
there are two horizons, however, Cor ley and Jacobson have shown that su- 
perluminal dispersion leads to an amplification of the particle production in 
the case of bosons. The analytic theory of this "black hole laser" process is 
quite complicated, so we provide some numerical results in the hope of aid- 
ing understanding of this interesting phenomenon. Specifically, we consider 
sonic horizons in a moving fluid. The theory of elementary excitations in 
a Bose- Einstein condensate provides an example of "super luminal" (Bogoli- 
ubov) dispersion, so we add Bogoliubov dispersion to Unruh's equation for 
sound in the fluid. A white-hole/black-hole horizon pair will then display 
black hole lasing. Numerical analysis of the evolution of a wave packet gives 
a clear picture of the amplification process. By utilizing the similarity of a 
radiating horizon to a parametric amplifier in quantum optics we also analyze 
the black hole laser as a quantum-optical network. 

1 Introduction 

One important contribution of the theoretical study of black hole analogues has been 
to help clarify the derivation of the Hawking effect [TJ EJ [3] . This in turn led to a 
study [1] of Hawking radiation in a more general context, one that involves, among 
other features, two horizons. The results of [I] are of particular interest because they 
offer a scenario, perhaps realizable in a black hole analogue, in which the Hawking 
radiation is amplified. Let us begin by describing the background to these ideas. 

There is an apparent contradiction in Hawking's semiclassical derivation of black 
hole evaporation [5] , in that the radiated fields undergo arbitrarily large blue-shifting 
in the calculation, thus acquiring arbitrarily large masses, which contravenes the 
underlying assumption that the gravitational effects of the quantum fields may be 
ignored. This is known as the trans-Plankian problem [HI [7J HI El E]. A similar 
issue arises in condensed matter analogues such as the sonic black hole [H [1] , and 
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Figure 1: Sonic horizons in a moving fluid, in which the speed of sound is 1. The 
velocity profile of the fluid, v(z), attains the value —1 at two values of z; these 
are horizons for sound waves that are right-moving with respect to the fluid. At 
the right-hand horizon right-moving waves are trapped, with waves just to the left 
of the horizon being swept into the supersonic flow region v < — 1; no sound can 
emerge from this region through the horizon, so it is reminiscent of a black hole. 
At the left-hand horizon right-moving waves become frozen and cannot enter the 
supersonic flow region; this is reminiscent of a white hole. 



to make our discussion concrete from the outset we shall consider sonic horizons in 
one-dimensional fluid flow (see Fig. [TJ. 

When the velocity profile of the fluid is as shown in Fig. [T] two horizons are 
formed for sound waves that propagate to the right with respect to the fluid. The 
horizon on the right of the supersonic flow region v < — 1 behaves like a black hole 
horizon for right-moving waves, while the horizon on the left of the supersonic flow 
region behaves like a white hole horizon for these waves. In [8] Unruh showed that 
in such a system, with some reasonable simplifying assumptions, the equation for a 
small perturbation cf) of the velocity potential is 

(d t + d z v)(d t + vd z )<j ) -d 2 z <j ) = 0. (1) 

In terms of a new coordinate r defined by 
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The results for quantum black holes El EH] apply equally well here, so it follows 
that each horizon will produce a thermal spectrum of phonons with a temperature 
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Figure 2: Hawking phonons in the fluid flow of Fig. [TJ Real phonons have positive 
frequency in the fluid-element frame and ^ shows that for right-moving phonons 
this frequency {u — vk) is = k. Thus in the subsonic-flow regions u (conserved 
for each ray) is positive, whereas in the supersonic-flow region it is negative; k is 
positive for all real phonons. The frequency in the fluid-element frame diverges at 
the horizons — the trans-Plankian problem. 



determined by the quantity that corresponds to the surface gravity at the horizon, 
namely the absolute value of the slope of the velocity profile: 
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The trajectories of the created phonons are readily drawn (see Fig. [2]), but 
they are formally deduced from the dispersion relation of the sound equation ([!]). 
Geometrical acoustics applied to ([T| gives the dispersion relation 



vk = ±k 



and the Hamilton equations 
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The left-hand side of ^ is the frequency in the frame co- moving with a fluid element, 
whereas u is the frequency in the laboratory frame; the latter is constant for a 
time- independent fluid flow ("time-independent Hamiltonian" du/dt = duj/dt = 0). 
Since the Hawking radiation is right-moving with respect to the fluid, we clearly 
must choose the positive sign in ^ and hence in Q also. By approximating v(z) as 
a linear function near the horizons we obtain from Q and ^ the ray trajectories 
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Figure 3: A "classical" picture of fluid atoms. Each atom is labeled and a is the 
equilibrium separation. 

of Fig. [2j The disturbing feature of the rays is the behavior of the wave vector k: 
at the horizons the radiation is exponentially blue-shiftecQ leading to a diverging 
frequency in the fluid-element frame. These runaway frequencies are unphysical 
since ([I]) asserts that sound in a fluid element obeys the ordinary wave equation 
at all wavelengths, in contradiction with the atomic nature of fluids (see Sec. [2]). 
Moreover the conclusion that this Hawking radiation is actually present in the fluid 
also assumes that ([T]) holds at all wavelengths, as exponential blue-shifting of wave 
packets at the horizon is a feature of the derivation [SHUdU]- Similarly, in the black- 
hole case the equation used in the calculation [9J [TO] does not hold at arbitrarily 
high frequencies because it ignores the gravity of the fields. For the black hole, 
a complete resolution of this difficulty will require knowledge of the gravitational 
physics of quantum fields, i.e. quantum gravity, but for the dumb hole the physics is 
available for a more realistic treatment. The issue to be addressed is the dispersion 
relation for sound at high frequencies p], and a consideration of this will lead us to 
the black hole laser. 



2 Dispersion 

In reality one would expect the dispersion relation for sound to differ from the linear 
formula ^ when the wavelength is of the order of the distance between the fluid 
atoms. What kind of dispersion relation should we expect at these wavelengths? 
Let us naively picture the fluid atoms as occupying well-defined positions with a 
separation distance a (see Fig. [3]). The equation for small oscillations in this chain 
of atoms is the wave equation with a discretized second-order spatial derivative: 

<9 t 2 0fc - — (4> k+1 - 20 fc + fc _i) = 0. (6) 
cr 

A partial differential equation that serves as an approximation to ^ is obtained by 
substituting the following relation, familiar from numerical analysis: 

\ - 2<j> k + fc _O = d 2 z 4> + ^k 4 + 0(a 4 ). 
er 12 

For small a we thus obtain the approximate dispersion relation 

^ - k 2 ~ £* 4 . (7) 



We trust the reader will not be too distracted by our mixing of acoustic and optical terminology. 
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As the last term in (|7j) is negative we find that the group velocity ^ of waves 
decreases with decreasing wavelength, i.e. we have subluminal or normal dispersion. 

A nonlinear dispersion relation to = uj(k) such as limits the amount of blue- 
shifting that occurs at a horizon, so the embarrassment of infinite frequencies disap- 
pears; there then remains the question of whether one can still derive the Hawking 
effect. Unruh pQ modified the sound equation ([I]) to incorporate subluminal disper- 
sion and showed numerically that the prediction of thermal radiation at temperature 
@ still holds good, regardless of the details of the high-frequency behavior. Indeed, 
analysis of generalized models shows that the precise form of the dispersion relation 
at high frequencies has no effect on the particle production by a single horizon, at 
least to lowest order [21 E] . Unsurprisingly, such a general conclusion involves some 
assumptions about the short- distance physics at the horizon, and since Planck-scale 
physics is unknown it remains possible that conditions at this scale conspire to 
invalidate the evaporation result for real black holes [3]. 

We now turn to the main purpose of this Section, namely the introduction of a 
dispersion relation that, together with the horizon arrangement of Fig. [TJ leads to 
black hole lasing. 



2.1 Bogoliubov dispersion 

The dispersion relation used by Corley and Jacobson in the original treatment of 
the black hole laser pE] is well known from the theory of elementary excitations 
in a Bose-Einstein condensate [TT4 [12] . In this formalism the bosonic atoms are 
described by an a field ip but most of the atoms are assumed to form a condensate 
with macroscopic wave function 

V>o = Vp~oe* So . (8) 

The deviations of ijj from the mean field ipo are given by a quantum field <fi of small 
fluctuations: 

It is then found that <j> satisfies the Bogoliubov-de Gennes equation 

ihd t <p = (- — V 2 - a + 2mc 2 ] } + mc 2 ft, (9) 
\ 2m J 

where [i is the chemical potential and c is the speed of sound. The field <p is given 
by the usual mode expansion 

= ^ {u v a u + v* u al) . (10) 

For plane- wave modes 

Uv = ue ik - r - lut , v v = ve^- r - iuJt , 
Q and its complex conjugate give 



vj \ -mc 2 + fi - 2mc z J \v 
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The condition for this matrix equation to have a solution is 



fh 2 k 2 2 \ / h 2 k 2 2 , \ 2 4 

u + 2mc — nw \ \ V U — 2mc" — nw + m c = 0, 

\ 2m J \ 2m ' ) 

and assuming the chemical potential \i is equal to the local energy of the condensate, 
mc 2 , this is the Bogoliubov dispersion relation 

<f = M + l?-k\ (11) 
4m 2 



From (11) we see that the group velocity increases with decreasing wavelength, 
so it is an example of superluminal or anomalous dispersion. Comparison with ^ 
shows that the purely quantum term in (11 ) corresponds to an imaginary separation 



distance between the atoms, a reflection of the inapplicability of the classical picture 
of Fig. [3] to a condensate. 

Let us now add Bogoliubov dispersion to our model system of Fig. [Tj The 
appropriate generalization of ([!]) is 

(d t + d z v)(d t + vd z )ct> - d 2 z <P + ^d 4 z <P = 0, (12) 

where k c is a constant; this give the dispersion relation 

I 

uo-vk = ±kJl + — , (13) 
so that we have Bogoliubov dispersion in a moving fluid. As with ^ we need only 



consider the positive sign in ( 13 ) since our interest is in phonons that propagate to 



the right with respect to the fluid. 

How does the presence of the /c c -dependent term in (13) alter the ray trajectories 
of Fig. [2]? We numerically solve the Hamilton equations resulting from ( 13 ) (with the 
positive sign) and obtain the trajectories shown in Fig. |4j The qualitative behavior 
of the rays in Fig. [4] is easily understood from the superluminal character of the 



dispersion relation (13). Consider the ray approaching the white hole horizon from 
the left: the ray moves at the group velocity and is blue-shifted as it nears the 
horizon; this increases the group velocity until the ray traverses the horizon into 
the supersonic-flow region. Similarly, the ray being swept towards the white hole 
horizon from the right is blue-shfted until the group velocity exceeds the speed of 
the flow, allowing the ray to move away from the horizon. The behavior at the black 
hole horizon is simply the time reverse of that at the white hole horizon. Thus, the 
ray in the supersonic-flow region bounces back and forth between the horizons while 
the ray incident from the left passes through the system. 

In Fig. [4] we have actually solved for the two rays using the same positive value of 
lu, the frequency in the laboratory frame, which is constant throughout the motion. 
It is then the case that the ray propagating through the system has positive frequency 
uj — v k in a frame co-moving with the fluid, whereas the bouncing ray has negative 
frequency in the fluid frame and does not therefore represent a real phonon. (A real 
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Figure 4: Numerical solution of ray trajectories for the dispersion relation ( 13 ). Both 
rays are right-moving with respect to the fluid and have the same positive value of 
uj. The ray that is confined to the supersonic-flow region has negative frequency in 
a frame co-moving with the fluid, whereas the ray that passes through the system 
from left to right has positive frequency in the fluid frame. 



phonon in the supersonic-flow region that is swept towards the white hole horizon 
has negative uj and exhibits the same bouncing behavior.) We have chosen the 
same value of uj for both rays in order to discuss particle creation. As the ray 
from the left reaches the horizon geometrical acoustics breaks down, allowing a 
wave packet centered on this ray to undergo mode conversion wherein it develops 
a component corresponding to another branch of the dispersion relation having the 
same conserved value of uj. Thus a wave packet centered on the incoming ray in Fig. [4] 
will develop a component between the horizons that is centered on the bouncing ray 
in Fig. [4], which has negative frequency in the fluid frame. In the quantized theory 
it is a positive frequency in the fluid frame that is associated with an annihilation 



operator through a mode expansion such as ( 10 ). That fact that a positive- frequency 
wave packet develops a negative-frequency component implies a mixing of creation 
and annihilation operators in the quantum theory and this means there is particle 
creation by the system [5j [9j [TUl [I] ■ 

The significance of the frequency in the fluid frame is apparent from the inner 
product 

hKdt + vd^-Mdt + vd^dz (14) 



which is conserved for solutions of (12), i.e. <9 t (0i,</> 2 ) = 0. The operator d t + vd z 
is the time derivative in the fluid frame and so is associated with the frequency in 
this frame. It is easily shown a wave packet made up of plane-wave solutions with 
positive (negative) frequency in the fluid frame has positive (negative) conserved 
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Figure 5: Schematic picture of the evolution of a positive-frequency wave packet <p- in . 
Lines that begin at the horizons denote wave packets created by mode conversion. 
All wave packets follow one of the ray trajectories of Fig. |4} The dotted lines indicate 
modes that are empty in this evolution but that feature in the connection formulae 
for modes at the horizons. 



norm (0, 0). Since it is the frequency in the fluid frame that is important for particle 
production and for the inner product (14), the terms positive/negative frequency 
will hereafter refer to the fluid frame. 

The above considerations lead to the schematic picture of Fig. [5] for the evolu- 
tion of an incident positive-frequency (and positive-norm) wave packet 0i n . At the 
white hole horizon the packet 0j n undergoes mode conversion, producing a negative- 
frequency component <j)^ R that propagates to the right together with the positive- 
frequency part (f>+. When these packets reach the black hole horizon there is further 
mode conversion and a positive- frequency packet 0^ exits the system leaving behind 
a red-shifted negative-frequency packet (jy_ L that propagates to the left towards the 
white hole horizon. On reaching the white hole horizon 0*^ generates through mode 
conversion a positive-frequency packet <j)\ , which joins the blue-shifted negative- 

(2) (2) 

frequency packet (fr_ R in propagating towards the black hole horizon. When 0+ and 

(2) 

(jr_ R reach the black hole horizon mode conversion occurs again and there is another 

(2) (2) 

output of a positive-frequency packet <% ut , leaving a negative-frequency packet <f>_ L 
between the horizons. The process continues as outlined. Now the total norm of 
this evolving solution is conserved, so in each instance of mode conversion amounts 
of positive- and negative-frequency wave packet are created that have equal and 
opposite norm. Thus with each bounce the norm of the negative-frequency packet 
between the horizons increases in magnitude; indeed, it increases exponentially since 
the norms of the wave packets created by mode conversion are in proportion to the 
norm of the incident packet. The size of the negative frequency packet generated by 
0i n determines the particle production so the number of particles produced increases 
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exponentially with time — this is the black hole laser [I]. 

To quantify the above description one can utilize the connection formulae ob- 
tained by Corley and Jacobson for modes at the horizons [I]. Let us assume that 
the slope of the velocity profile is —a at the white hole horizon and 7 at the black 
hole horizon. Then the relation between modes with the same conserved value of u> 
at the white hole horizon is (see Fig. [5]): 



a;— Lin 



[2 sinh {iru/a)]' 
- [2 sinh (iru/a)]' 



1/2 ( Mj/ZaAl) 



+ e 



-R ' 



-1/2 ^-7ra;/2a^(l) + ^/Za^l) 



R ' 



(15) 
(16) 



with similar formulae for the other modes that meet at this horizon in Fig. |5j At 
the black hole horizon the connection formulae for modes is (see Fig. [5]): 



</fi ut = [2 sinh (W7)r 1/2 (V^tfj + e™l**4>£ 



3-L ) ' 



- [2 sinh (™/7)r 1/2 i e-^^l + e^ 2 ^ (1) 



j uj-L 



ui—L ; > 



(17) 
(18) 



with similar formulae applying for the other modes that interact at this horizon 
in Fig. 5 For modes that travel between the horizons <P^-r an d <P^1l) we 

must in addition specify how their values at the two horizons are related to each 
other. This can be done using a WKB approximation to propagate the modes from 
one horizon to the other; the result will be a simple phase relationship between the 
values at each horizon: 
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(19) 
(20) 
(21) 



where BH (WH) refers to the black hole (white hole) horizon. Using the above 
formulae as building blocks one can calculate the evolution of a wave packet and 
hence compute the number of created particles [1] ; in Sec. [4] we shall address this 
task in a different guise. First, a numerical solution for an evolving wave packet will 
provide an enlightening picture of the process outlined in Fig. [5] 



3 Numerical results 

Let us summarize what we expect to see if we propagate a positive-frequency wave 
packet towards the white hole horizon. According to Fig. [5] this positive-norm wave 
packet will plough through the system and emerge out of the black hole horizon 
(0out) somewhat amplified, leaving behind a small negative-norm packet (<fi_ L ) such 
that the total norm is conserved. After a period of time equal to that required 
by the negative-frequency ray in Fig. [4] to propagate from the black hole horizon 
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Figure 6: Numerical solution for a positive-frequency wave packet that has entered 
the horizon system from the left. The packet is centered on the positive- frequency 
ray in Fig. [4] so it has propagated through the system (0^) and begun to emerge 
on the left-hand side due to periodic boundary conditions. A sequence of positive- 
frequency wave packets are seen to follow 0ouL i n nne with expectations based on 
Fig. [5} 



to the white hole horizon and back again we expect another, smaller, positive-norm 
packet (0o Ut ) to emerge from the system. At regular intervals there should be further 
outputs of positive-norm packets (0o Ut , etc.) and a build-up of the negative-norm 
packet between the horizons. The increasing amplitude of the negative-norm packet 
will cause an increase in the size of the positive-norm outputs produced by mode 
conversion, and this increase should in fact be exponential. 

In Fig. [6] we show the result of a numerical evolution of a wave packet centered 
on the positive-frequency ray in Fig. |4j Periodic boundary conditions are in force 
in Fig. 6] and the incident packet has ploughed through the system (4>^t) an d is 
circling round on the left-hand side where it is about to hit the white hole horizon a 
second time. The successive outputs of positive-norm packets (<%ut> e tc) are clearly 
visible; five separate emissions are distinguishable, corresponding to the negative- 
norm packet bouncing off the black hole horizon five times. As the total amount of 
positive-norm packet increases with each output it is balanced by a growing amount 
of negative-norm packet between the horizons, leaving the total norm unchanged. 

The distance between the horizons has been set quite small relative to the length 
of the incident wave packet in order to increase the number of bounces that occur 
before <f>^ t collides with the white hole horizon. This makes it difficult to see clearly 
the motion of the negative-frequency wave packet that is confined to the supersonic- 
flow region. Nevertheless alternate red- and blue-shifting is observed between the 
horizons as the solution evolves and the time between successive positive-norm out- 
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puts is equal to the time between visits to the black hole horizon of the ray on which 
the negative-frequency packet is centered (see Fig. [4]). 

Although the last output in Fig. |(3| seems to be somewhat larger than 
there is no clear sign of an exponential increase in the size of the output packets 
after this number of emissions. The analytic formulae for the evolution show that 
about ten outputs would be required in this simulation before there is a significant 
increase in their size, as we shall see in the next Section. 

Before we turn to an analytic treatment of the black hole laser and the calcu- 
lation of the particle production, we briefly mention the equally remarkable results 
when the field is fermionic rather than bosonic The conserved norm for a 
fermionic field is positive definite, so both positive- and negative-frequency packets 
have positive norm. From this fact alone follows the evolution of a positive-frequency 
fermionic wave packet. When the initial wave packet ploughs through the system 
as in Fig. [6] it is diminished in amplitude in order to conserve the total norm, since 
the negative-frequency packet left behind also has positive norm. As subsequent 
positive-frequency packets are emitted conservation of norm requires the amount of 
negative frequency packet between the horizons to decrease rather than increase as 
in the bosonic case. In the limit of infinite time the amount of negative frequency 
packet goes to zero and therefore so does the number of particles produced — there 
is complete suppression of the Hawking effect. This can be understood as a filling 
up of all the allowed fermionic states between the horizons by the Hawking particles 
created in this region jl]. 



4 Black hole amplifier 



The Bogoliubov transformation that underlies particle creation by a black hole is 
similar to that which describes the action of a parametric amplifier or a phase- 
conjugating mirror in quantum optics [131 EH]- An optical device will transform a set 
of incident quantum light modes, described by annihilation operators {eii, . . . , a n } 



and their Hermitian conjugates {a[, . 
operators {a[, . . . , a' n } and {a'/, . . . , a' n 
is given by the scattering matrix S: 

/a'A 



/t 



. .,at}, into outgoing modes, described by 
'}. If the device is linear this transformation 



S 



(22) 



The S'-matrix is constrained by the requirement that the creation and annihilation 
operators satisfy the usual commutation relations; this imposes the quasi-unitarity 
condition 

SGS^ = G, G = ( T Q _° 7 1 . (23) 
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Matrices S satisfying (23 ) form the group U (n, n) 
is 

/ cosh £ 



An example of such a matrix 



S 






V sinh £ 





cosh£ 
sinh£ 






sinh£ 
cosh£ 




sinh £\ 



cosh£ / 



(24) 



which is the scattering matrix for a parametric amplifier or a phase-conjugating 
mirror [T3J HI] . In this case there are two input and output modes and the device 
creates or annihilates pairs of photons. The energy for the pair production or the 
reservoir for pair annihilation is provided by the pump process of the amplifier, 
which determines the parameter £. For a black hole we may define £ by 



tanh£ 



3 — 



(25) 



where 7 is the surface gravity at the horizon and u is the frequency of the quantum 



modes under consideration. Then (24) gives the Bogoliubov transformation for the 



Hawking effect [U [TU] , which we write concisely as 



/ cosh £ sinh £ 
y sinh £ cosh £ 



a 1 



(26) 



For a moving fluid, (25) and (26) describe phonon creation by a black hole 



horizon, where 7 denotes the slope of the velocity profile v(z) at the horizon. A 
white hole horizon in the fluid, where the slope of v(z) is —a, is the time reverse of 
a black hole horizon. The "surface gravity" is a and the roles of input and output 
modes are reversed. Thus, defining a parameter £ by 



tanh ( 



-ITUj/oi 



the white hole horizon generates the Bogoliubov transformation 



a 1 



( cosh £ — sinh £ 
1— sinh£ cosh^ 



(27) 



(2f 



In the black hole laser the particle production by the pair of horizons is not given 



by (26) and (28) because the horizons interact due to the superluminal dispersion. 



We can represent the black hole laser as a quantum-optical device by expressing the 
amplification process of Fig. [5] in terms of the annihilation operators of the modes, 
and using a more convenient notation (see Fig. [7]). In the nth amplification step we 
have a pair of input modes a +n and a_ n at the white hole horizon that generates an 
output mode a' +n on the right-hand side of the black hole laser as well as a mode 
a'_ n that is trapped between the horizons where it serves as a further input to the 
device, i.e. 

a'_ n = a-(„+i). (29) 
The picture is that of a network of amplifiers, as outlined in Fig. [7} 
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Figure 7: The black hole laser as a network of amplifiers. The mode operators 
with negative subscripts correspond to the Hawking radiation trapped between the 
horizons, while those with positive subscripts describe the Hawking radiation present 
outside the horizons. An input consisting of a mode on the left of the black hole 
laser and a mode between the horizons produces an output mode on the right of 
the black hole laser and another output mode between the horizons, as in Fig. |5j 
The process repeats, with the output mode between the horizons acting as an input 
mode for the next cycle, i.e. d'_ n = d_( n+ i). 



The output of the black hole laser is completely specified by the matrix B that 
effects the Bogoliubov transformation for each step in Fig. [7J 

B ( *r) ■ (30) 



-ti 



There are four distinct processes incorporated in the matrix B, allowing it to be 
expressed as 

B = B A BzB 2 B x . (31) 



First there is the pair production by the white hole horizon, described by (28); we 
therefore have 

Bi=( C ° Sh is " Sin A C V (32) 



V — sinh ( cosh ( 

Next, the two modes produced by the white hole horizon propagate to the black 
hole horizon, as in Fig. |5j this imparts separate phases and 9 + to each mode 
operator, so we write 



B2=[ e * + ). (33) 



At the black hole horizon the modes are transformed according to (26), giving 

Bs=( C ° Sh h l Si tiY (34) 
V sinh 4 cosh £ / v ' 
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Finally the mode trapped between the the horizons propagates to the white hole 
horizon, acquiring a phase 9q\ 



Ba 



e -i0 O o 
1 



(35) 



The particle content after m amplification steps is now easily calculated from ( 30 ) 
and (29 ), using the value (31 )-(35 ) for B. It is worthwhile to consider the restrictions 
on the 2x2 matrix B resulting from the fact that it determines a U(2, 2) scattering 
matrix S (see (30) and (22)). Constructing S from B and imposing (23) one finds 
that the components of B must satisfy 



\Bn 



\B 



12 



1. 



\B 



21 



IB- 



22 



1, B* n B 21 - B* 12 B 22 = 0, (36) 



which is the statement that B is a U(l, 1) matrix. The constraints (36) impose the 
following structure on B: 



B 



* , i> real, 



\v\ 



1. 



(37) 
(38) 



between the horizons and a+ m for the mode outside the horizons are given by 



Using (|30|), (|29|) and (|37|) we find after m cycles that d'_ m for the mode trapped 

(39) 



+ v*J2e i( - m - l+1 ^fji m - l a< +l , 



i=i 



m— 1 



e im ^vn m ~ x a^ + \v\ 



e l{m - l+ ^^ m - l - l a ] +l + e l ^n*ai m . 



(40) 



l=i 



The condition that we have vacuum before the formation of the black hole laser 
gives 

6_i|0) = 0, a+;|0> = 0, / = l,...m. 

We may write the particle content for the two output modes in terms of |/i| 2 because 
of (38). The particle number for the trapped mode after m cycles is 

(AL m ) = (a'l m a'_ m ) 

1 2m y 



1 2 



2(m-l) 



I Ml 



I/' 



i=i 

2m 



l 



(41) 



whereas the number of particles outside the horizons is 

(N +m ) = (a'l m a +m ) 

= |z/| 2 |/i| 2(m_1) 

= H 2m (i-H- 2 )- 
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(42) 



Figure 8: The number of Hawking particles outside the horizons versus the number 



of amplification cycles, as given by formula (42). 



The explicit form of |/i| 2 is found from (31)-(35) and (37) to be 



iMl 



1 + cosh(2£) cosh(2C) - cos(#+ + 6 J) sinh(2£) sinh(2C) 



(43) 



Note that the phase 9q does not appear in (43) and thus has no effect on the particle 



content; this is because 9q is acquired by a mode traveling to the white hole horizon 
where it interacts with a mode from outside the device, and the latter has no phase 
sensitivity. The combination 9 + + 9_, however, is the phase difference between the 
two modes that propagate across the device from the white hole horizon and this 



affects how these modes interact at the black hole horizon. We see from (41)-(43) 
that for frequencies such that cos(# + + 0-) = 1 the particle production is maximally 
suppressed; indeed if the horizons are symmetric (£ = () there is no production at 
this frequency. On the other hand maximal amplification of the Hawking radiation 
occurs at frequencies for which cos(# + + 6 J) = 0, giving a WKB-type condition 



9+ + 9- = 2tt in + 



n an integer. 



In Fig. 8 we use formula (42) to plot the particle number outside the horizons as 



a function of the number of amplification cycles. The slopes of v(z) at the horizons 
are those used for the numerical example in Sec. [3] and the frequency chosen is that 
on which the incident pulse in Fig [6] is centered. A WKB approximation is used to 
estimate the phases 9 + and 6L. Fig. [8] shows the consequence, in terms of Hawking 
radiation, of the wave-packet evolution in Fig. [6| As discussed in Sec. [3j the stream 
of positive-frequency pulses emerging from the black hole horizon in Fig. [6] implies 
a growth in particle production, as in Fig. [8j The solution domain in Fig. [6] is 
only large enough to contain the result of five amplification cycles. In this time 
the positive-frequency pulses trailing are of approximately the same size; this 
implies an approximately linear increase in the particle number, as is the case in 
Fig. [8] with n < 5. In order for the exponential nature of the increase in particle 
number to be clearly seen in a wave-packet evolution, one would require at least ten 
positive-frequency outputs in the evolution, as one sees from Fig. [8] 

In summary, the black hole laser is a dramatic quantum-vacuum effect and pro- 
vides an interesting example of the interplay of Hawking radiation and nonlinear 
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dispersion. Any cosmological significance of this process remains highly specula- 
tive jl], but in contrast a laboratory example of black hole lasing is not nearly so 
far-fetched. Bose-Einstein condensates have the properties required for the con- 
struction of a sonic black hole laser, however formidable the practical difficulties 
may be, and this is not the only current possibility |16j . As the study of black hole 
analogues continues the prospect of exponentially amplified Hawking radiation is 
tantalizing for both theorists and experimentalists. 
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